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ABSTRACT 

The Fermi-LdLYgo, Area Telescope (EAT) First Source Catalog (IFGL) was released in February 2010 and the 
Fermi-hAI 2-Year Source Catalog (2FGL) appeared in April 2012, based on data from 24 months of operation. 
Since their releases, many follow up observations of unidentified gamma-ray sources (UGSs) were performed 
and new procedures to associate gamma-ray sources with potential counterparts at other wavelengths were 
developed. Here we review and characterize all the associations as published in the IFGL and 2FGL catalog 
on the basis of multifrequency archival observations. In particular we located 177 spectra for the low-energy 
counterparts that were not listed in the previous Fermi catalogs, and in addition we present new spectroscopic 
observations of 8 y-ray blazar candidates. Based on our investigations, we introduce a new counterpart category 
of “candidate associations” and propose a refined classification for the candidate low-energy counterparts of the 
Fermi sources. We compare the IFGL-assigned counterparts with those listed in the 2FGL to determine which 
unassociated sources became associated in later releases of the Fermi catalogs. We also search for potential 
counterparts to all the remaining unassociated Ei^rm/sources. Finally, we prepare a refined and merged list of 
all the associations of the IFGL plus 2FGL catalogs that includes 2219 unique Fermi objects. This is the most 
comprehensive and systematic study of all the associations collected for the y-ray sources available to date. 
We conclude that 80% of the Fermi sources have at least one known plausible gamma-ray emitter within their 
positional uncertainty regions. 

Subject headings: methods: statistical - galaxies: active - quasars: general - surveys 


1. INTRODUCTION 

Despite the large progress in y-ray source localiza¬ 
tion made by the Fermi Large Area Telescope (EAT) 
(Atwood et al. 2009) coupled with the improvements in our 
knowledge of the diffuse Galactic y-ray emission (e.g., 
Moskalenko et al. 2007; Abdo et al. 2009), the positional un¬ 
certainties of the Fermi sources are still large with respect 
to typical localization precisions at X-ray and lower ener¬ 
gies, making a significant fraction of the sources in the y- 
ray sky yet unknown (Abdo et al. 2010a; Nolan et al. 2012). 
About 43% of the Fermi sources detected in the Fermi- 
Large Area Telescope (EAT) First Source Catalog (IFGL 
Abdo et al. 2010a) were listed as UGSs while there were 
-33% y-ray objects unassociated in the Fermi-LKF 2-Year 
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Source Catalog (2FGL Nolan et al. 2012). 

To decrease the number of unidentified gamma-ray sources 
(UGSs) many methods based on multifrequency approaches 
or statistical analyses have been recently adopted. Radio 
follow up observations of the Fermi UGSs have already 
been performed (e.g., Kovalev 2009; Hovatta et al. 2012; 
Petrov et al. 2013; Hovatta et al. 2014; Schinzel et al. 2014) 
and the Swift X-ray survey for all the UGSs listed 
in the Fermi source catalogs is still on going^^ 
(e.g., Mirabal & Halpern 2009; Paggi et al. 2013; 
Takeuchi et al. 2013; Stroh et al. 2013; Acero et al. 2013). 
Additional X-ray observations performed with Chandra 
and Suzaku improved our knowledge on the UGSs (e.g., 
Maeda et al. 2011; Cheung et al. 2012; Kataoka et al. 2012; 
Takahashi et al. 2012). Statistical studies based on 
the y-ray source properties also allowed us to recog¬ 
nize the nature of the potential counterparts for the 
UGSs (e.g., Ackermann et al. 2012; Mirabal et al. 2012; 
Hassan et al. 2013; Doert & Errando 2014). More¬ 

over a tight connection between the infrared (IR) sky 
seen by Wide-Field Infrared Survey Explorer (WISE 
Wright et al. 2010) and the Fermi one has been re¬ 
cently discovered for blazars, the rarest class of active 
galaxies (Massaro et al. 2011a; D’Abrusco et al. 2012; 
D’Abrusco et al. 2013). These works greatly decreased 
the fraction of UGSs with no assigned counterpart at 
low energies (Massaro et al. 2012a; Massaro et al. 2012b; 
Massaro et al. 2013a). More recently also low frequency 
radio observations (i.e., below -1 GHz) revealed a new 
spectral behavior that allowed us to search for blazar-like 
sources lying within the positional uncertainty regions of the 
UGSs (Massaro et al. 2013b; Nori et al. 2014). Additional 

http://www.swift.psu.edu/unassociated/ 



2 


F. Massaro et al. 2014 


studies have been also carried out with near-infrared obser¬ 
vations (Raitieri et al. 2014) as well as in the sub millimeter 
range (e.g., Giommi et al. 2012; Lopez-Caniego et al. 2013). 
Optical spectroscopic campaigns also have been cru¬ 
cial to disentangle and/or confirm the natures of 
the low-energy counterparts selected with different 
methods (e.g., Masetti et al. 2013; Shaw et al. 2013a; 
Shaw et al. 2013b; Paggi et al. 2014; Massaro et al. 2014b; 
Landoni et al. 2014). 

The main aim of the analysis presented here is to con¬ 
firm the current status of the associations for the IFGL and 
the 2FGL y-ray sources. Thus we prepare a single list, cre¬ 
ated by merging the unique sources in the IFGL and in 
the 2FGL catalogs. On the basis of the information re¬ 
ported in the First Fermi-LNI Catalog of Sources Above 
10 GeV (IFHL Ackermann et al. 2013) and in both the First 
and the Second Catalog of Active Galactic Nuclei Detected 
by Fermi (ILAC and 2LAC, respectively Abdo et al. 2010b; 
Ackermann et al. 2011a), together with an extensive litera¬ 
ture search in multifrequency archives, we verify the cur¬ 
rent status of the y-ray associations. In addition, we also 
present new optical spectroscopic observations of eight y-ray 
blazar candidates selected during our multifrequency study. 
The result is a refined and merged list of all the associations 
for the lFGL-plus-2FGL catalog with updated information 
on the low-energy counterparts including specific multifre¬ 
quency notes. 

The paper is organized as follows: in Section 2 we dis¬ 
cuss the Fermi procedures to assign identifications and asso¬ 
ciations with potential counterparts at low energies, and the 
categories of y-ray source associations. Section 3 describes 
the content and the main properties of both the IFGL and 
the 2FGL catalogs. In Section 4 we introduce an updated 
y-ray source classification for the counterparts of the Fermi 
detected objects. Details on the correlation with multifre¬ 
quency databases and catalogs are in Section 5. In Section 6 
we present the refined version of the associations listed in the 
Fermi catalogs (hereinafter designated IFGLR and 2FGLR, 
respectively, for brevity). The analysis of the infrared colors 
for the sources associated in the IFGL is given in Section 7 
and new optical spectroscopic observations of y-ray blazar 
candidates are presented in Section 8. In Section 9 we com¬ 
pare our results with those achieved from statistical analyses 
performed on the UGSs listed in both the IFGL and the 2FGL. 
In Section 10 we speculate on both the radio-y-ray and the 
IR-y-ray connections for the Fermi blazars to check the con¬ 
sistency of sources classified as blazar candidates with these 
multifrequency behaviors. The summary and conclusions are 
given in Section 11. We use cgs units unless stated otherwise. 
Spectral indices, a, are defined by flux density, Sy oc and 
WISE magnitudes at [3.4], [4.6], [12], [22] jum (i.e., the nom¬ 
inal bands) are in the Vega system. For numerical results a 
flat cosmology was assumed with Hq = 67.3 km s“^ Mpc“\ 
Qm = 0.315 and Qa = 0.685 (Planck collaboration 2014). 
The astronomical survey and source class acronyms used 
in the paper are listed in Table 1. 

2. CATEGORIES OF GAMMA-RAY SOURCE ASSOCIATIONS 

In all the Fermi catalogs there is an important distinc¬ 
tion between identification of low-energy counterparts for 
the Fermi sources and association. Identification is based 
on i) spin or orbital periodicity (e.g., pulsars, binary 
systems) or on ii) correlated variability at other wave¬ 
lengths (e.g., blazars, active galaxies) or on iii) the con¬ 


TABLEI 

List of acronyms used in the paper. 


Acronym 

Term 

IFGL 

Fermi-LdsgQ Area Telescope First Source Catalog 

IFGLR 

IFGL Refined association catalog 

IFHL 

First Fermi-Large Area Telescope Catalog of Sources Above 10 GeV 

ILAC 

First Catalog of Active Galactic Nuclei Detected by Fermi 

IXSPS 

Deep Swift X-Ray Telescope Point Source Catalog 

2FGL 

Fermi-Large Area Telescope 2-Year Source Catalog 

2FGLR 

2FGL Refined association catalog 

2LAC 

Second Catalog of Active Galactic Nuclei Detected by Fermi 

2MASS 

Two Micron All-Sky Survey 

3C 

Third Cambridge Catalogue of Radio Sources 

6dFGS 

Six-degree-Field Galaxy Redshift Survey 

AGN 

Active Galactic Nucleus 

AGU 

AGN of Uncertain type 

AT20G 

Australia Telescope 20 GHz Survey 

ATNF 

Australia Telescope National Facility 

BZB 

BL Lac object 

BZQ 

flat spectrum radio quasar 

CRATES 

Combined Radio All-Sky Targeted Eight-GHz Survey radio catalog 

CSC 

Chandra Source Catalog 

CT 

Classification Tree method 

FIRST 

VLA Faint Images of the Radio Sky at Twenty-Centimeters 

GB6 

Green Bank 6-cm Radio Source Catalog 

GLC 

Globular Cluster 

HMB 

High Mass X-ray Binary 

KDE 

Kernel Density Estimation 

IRAF 

Image Reduction and Analysis Facility 

IRAS 

InfraRed Astronomical Satellite 

LAT 

Large Area Telescope 

LBA 

Australian Long Baseline Array 

LCSI 

LBA Calibrator Survey 

LORCAT 

Low-frequency Radio Catalog of flat-spectrum Sources 

LR 

Logistic Regression method 

MGPS 

Molonglo Galactic Plane Survey 

MSP 

Millisecond Pulsar 

NED 

NASA/IPAC Extragalactic Database 

NOV 

Nova 

NRAO 

National Radio Astronomy Observatory 

NVSS 

NRAO VLA Sky Survey Catalog 

OAGH 

Observatorio Astrofisico Guillermo Haro 

OAN 

Observatorio Astronomico Nacional 

PKS 

Parkes Southern Radio Source catalog 

PMN 

Parkes-MIT-NRAO Surveys 

PSR 

Pulsar 

PWN 

Pulsar wind nebula 

RBSC 

ROSAT Bright Source Catalog 

RFSC 

ROSAT Faint Source Catalog 

SDSS 

Sloan Digitized Sky Survey 

SED 

Spectral Energy Distribution 

SFR 

Star Forming - HII regions 

SNR 

Supernova Remnant 

SUMSS 

Sydney University Molonglo Sky Survey 

TEXAS 

Texas Survey of Radio Sources at 365 MHz 

UGS 

Unidentified Gamma-ray Source 

VLA 

Very Large Array 

VLSS 

VLA Low-Frequency Sky Survey Discrete Source Catalog 

USNO 

United States Naval Observatory 

XMMSL 

XMM-Newton Slew Survey 

WENSS 

Westerbork Northern Sky Survey 

WISE 

Wide-Field Infrared Survey Explorer 

WISH 

Westerbork in the Southern Hemisphere Source Catalog 

WSRT 

Westerbork Synthesis Radio Telescope 

WSRTGP 

WSRT Galactic Plane Compact 327-MHz Source Catalog 


sistency between the measured angular sizes in y-ray and 
at lower energies (e.g., supernova remnants). On the other 
hand, the association designation indeed depends on the re¬ 
sults of different procedures adopted in both Fermi catalogs 
(Abdo et al. 2010a; Nolan et al. 2012). These procedures are: 

1. The Bayesian Association Method: Initially applied 
to associate EGRET sources with flat-spectrum radio 
sources (e.g., Mattox et al. 1997; Mattox et al. 2001; 
Abdo et al. 2010a), this method assesses the probabil¬ 
ity of association between a y-ray source and a can¬ 
didate counterpart taking into account their local den¬ 
sities. Local density is estimated simply by counting 
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candidates in a nearby region of the sky. 

2. The Likelihood Ratio method: Used to search for possi¬ 
ble counterparts in uniform surveys in the radio and in 
X-ray bands, this procedure was originally proposed by 
Richter (1975) and subsequently applied by and modi¬ 
fied by de Ruiter, Willis & Arp (1977), Prestage & Pea¬ 
cock (1983), Wolstencroft et al. (1986) and by Suther¬ 
land & Saunders (1992). 

3. The logN-logS association method: This is a modi¬ 
fied version of the Bayesian method for blazars taking 
into account their log A-log^ (see Abdo et al. 2010b; 
Ackermann et al. 2011a, for more details). 

There are several differences between the source asso¬ 
ciations presented between the IFGL and the 2FGL cat¬ 
alogs mostly related to the y-ray analysis and to the im¬ 
provements achieved in the models of Galactic y-ray dif¬ 
fuse emission as well as in evaluation of the LAT response 
functions. Improving the Fermi source localization not only 
makes easier the search for low-energy counterparts but al¬ 
lows more accurate estimates of the association probability 
since they depend on the y-ray positional uncertainty (see 
e.g., Abdo et al. 2010b; Ackermann et al. 2011a). These dif¬ 
ferences also relate to updates of comparison surveys and cat¬ 
alogs used to search for the low-energy counterparts and on 
multifrequency observations performed on UGS samples. 

The use of catalogs/surveys that were previously not con¬ 
sidered for the Fermi catalog preparation, as for example 
the infrared all-sky survey performed by WISE and the low- 
frequency radio catalogs, revealing new connections between 
the low and the high-energy skies (e.g., Massaro et al. 2011a; 
D’Abrusco et al. 2013; Massaro et al. 2013b) has also de¬ 
creased the unknown fraction of the y-ray sky (e.g., 
Massaro et al. 2012b; Massaro et al. 2013a). This affects the 
Fermi associations because the methods adopted to assign 
low-energy counterparts depend on the source densities of the 
catalogs (Nolan et al. 2012; Ackermann et al. 2011a). A bet¬ 
ter estimate of the counterpart density leads to a better esti¬ 
mate of the false positive associations and of the association 
probability and this implies that a previously unassociated 
source could have an assigned counterpart in a new release 
of the Fermi catalog. 

Motivated by the changes occurring between the IFGL and 
the 2FGL associations here we propose to introduce, together 
with identified and associated sources, the category of ‘'can¬ 
didate associations”. These are y-ray sources that have a 
potential low-energy counterpart in a specific class of well- 
known y-ray emitters lying within the Fermi positional uncer¬ 
tainty region at 95% level of confidence and/or having angu¬ 
lar separations between the Fermi and the counterpart position 
smaller than the maximum one for all the associated sources 
of the same class. Candidate associations have the potential 
to be promoted to associations in the future (see Section 4). 

We considered as counterpart classes for the candidate as¬ 
sociations only a few kinds of sources: i) blazar-like sources, 
as defined in the following, ii) known pulsars (PSRs), iii) pul¬ 
sar wind nebulae (PWNe), supernova remnants (SNRs) and 
star forming and HII regions (hereinafter simply indicated as 
SFRs). The first two source classes constitute the two largest 
populations of known y-ray emitters while the other choices 
were considered because 15 out of the 44 SNRs and 10 out 
of the 63 PSRs and PWNe associated in the IFGL lie within 
SFRs and a similar situation occurs for the 2FGL. 


Our choice concerning the SFRs, especially for the most 
massive cases, is indeed motivated by the possibility that they 
could host type O and/or B stars, potentially associated with 
open star clusters, and this being relatively likely to have 
some members that have progressed to become a PSR, PWN, 
and/or SNR, all known classes of y-ray sources (see also 
Montmerle 1979; Montmerle 2009). 

In addition, considering SFRs as candidate associations is 
also motivated by the recent idea, partially supported by the 
Fermi discovery of y-ray emission arising from Eta Cari- 
nae (Abdo et al. 2010c), that shells detected in the radio, in¬ 
frared, and X-rays could mark termination shocks of stel¬ 
lar winds. These winds interacting with the stars’ natal 
molecular clouds can be considered as plausible sites to ac¬ 
celerate particles that emit y-rays (e.g., Araudo et al. 2007; 
Bosch-Ramon et al. 201; Rowell et al. 2010). Bowshocks are 
associated with many source classes such as PSRs, cata¬ 
clysmic variables, colliding wind binaries, cometary H II re¬ 
gions and particles accelerated therein can radiate up to the 
MeV-GeV energy range (see e.g. Del Valle & Romero 2012). 
However we remark that evidence of y-ray emission from 
colliding wind binaries associated with Wolf-Rayet stars was 
not found by Fermi, for a selected sample of 7 sources 
(Werner et al. 2013); although, as highlighted by the same au¬ 
thors, searching for such y-ray detections is extremely com¬ 
plicated due to the uncertainties of the diffuse emission in the 
Galactic plane. 

Some candidate associations lack multifrequency infor¬ 
mation as well as correct spectroscopic identifications but 
exhibit some characteristic features, typical of known y- 
ray emitters, mainly blazars and PSRs. They were not 
associated in the previous releases of the Fermi catalogs 
mostly because they were not listed in one of the appro¬ 
priate comparison catalogs of potential counterparts used 
to associated Fermi sources. However once follow up 
observations were performed (see e.g. Paggi et al. 2013; 
Masetti et al. 2013; Petrov et al. 2013; Massaro et al. 2014b; 
Landoni et al. 2014, for blazar-like candidates), then they 
could be promoted to associated sources in new releases of 
the Fermi catalogs. In addition, within this category we also 
considered Fermi sources having SFRs within their y-ray po¬ 
sitional uncertainty regions for which multifrequency obser¬ 
vations could reveal the presence of unknown PSRs, PWNe 
and SNRs, already counterparts of ~5% of the y-ray sources. 

We emphasize three advantages obtained from introducing 
the concept of candidate associations: 

1. Designating candidate association category allows us to 
flag potential counterparts that might be associated or 
identified in future releases of the Fermi catalogs, as 
previously occurred to 171 IFGL sources then recog¬ 
nized, classified and associated in the 2FGL. This also 
facilitates planning follow up multifrequency observa¬ 
tions, in particular spectroscopic ones, necessary to de¬ 
termine and/or confirm their nature. 

2. The category of candidate associations permits 
us to establish which is the number of Fermi 
sources not yet having potential counterparts. 
The precise knowledge of this number is ex¬ 
tremely relevant to set better constraints on 
dark matter scenarios (e.g., Zechlin et al. 2012; 
Mirabal et al. 2012; Mirabal 2013a; Mirabal 2013b; 
Berlin & Hooper 2014). 


4 


F. Massaro et al. 2014 


3. For all the methods adopted to associate sources in the 
Fermi catalogs, the resulting associations depend on the 
comparison catalog/survey used. The same low-energy 
counterpart may be associated with a Fermi source us¬ 
ing a certain catalog but not with a different survey, 
even when the positions in the two catalogs are the 
same. This occurs because the source density of the 
comparison catalogs, which is used to determine the 
probability of false positive associations, affects the as¬ 
sociation probability. Thus observing these candidate 
associations with follow up campaigns to determine 
their natures will permit us to add blazars and pulsars 
to the correct comparison catalogs and associate them 
in future releases of the Fermi surveys. 

To reiterate the importance of having candidate associa¬ 
tions we highlight the recent work of Schinzel et al. (2014) 
on radio very-long-baseline interferometry (VLBI) observa¬ 
tions of Fermi unassociated sources in the 2FGL (see also 
Petrov et al. 2013). In the final list of high-confidence as¬ 
sociations presented in Schinzel et al. (2014) more than 
90% of the sources are also detected in the major radio 
surveys such as the National Radio Astronomy Observa¬ 
tory (NRAO) Very Large Array (VLA) Sky Survey Cat¬ 
alog (NVSS; Condon et al. 1998) and the Sydney Univer¬ 
sity Molonglo Sky Survey (SUMSS; Mauch et al. 2003). 
These counterparts were not associated using the methods 
adopted in Fermi catalogs due to the high density of sources 
at low flux level in the above surveys. However thanks to 
the high resolution radio follow up observations, which rev¬ 
eled the presence of radio compact cores, it has been pos¬ 
sible to use a different selected sample of potential counter¬ 
parts to find new associations via the likelihood procedure 
(Schinzel et al. 2014; Petrov et al. 2013). Thus it is extremely 
important to have targets for follow up observations that, once 
performed, will allow us to insert the potential counterparts in 
the correct catalogs used to search y-ray associations. 

Identifications, associations and candidate associations’ of 
y-ray sources are labelled in both the refined version of the 
Fermi catalogs as the category of the y-ray source associa¬ 
tion. These are defined as: identified (I), associated (A), can¬ 
didate association (C) in the tables while sources for which 
no potential counterpart was found are assigned to the UGS 
category (U). 

3. THE Fermi CATALOGS 

3.1. The Fermi-Large Area Telescope First Source Catalog 

The IFGL catalog^^ contains 1451 sources detected and 
characterized in the 100 MeV to 100 GeV energy range. For 
each Fermi source listed therein, positional uncertainty re¬ 
gions at 68% and 95% level of confidence, results from spec¬ 
tral fits as well as flux measurements in five energy bands are 
reported (Abdo et al. 2010a). Firm identification or plausible 
association with a low-energy counterpart was provided for 
each IFGL source on the basis of the comparison with other 
astronomical catalogs and on the basis of the association pro¬ 
cedures previously mentioned (see Section 2 and Abdo et al. 
2010a for details). 

The published IFGL catalog lists several classes of coun¬ 
terparts: 63 pulsars (PSRs), 7 associated plus 56 identified 
(24 of them discovered thanks to the Fermi-LP^ observa¬ 
tions), 2 pulsars wind nebulae (PWNe) 44 y-ray sources cor- 

^^ http://fermi.gsfc.nasa.gov/ssc/data/access/lat/1 yr.catalog/ 


responding to 38 SNRs, 41 associated and 3 identified, 8 glob¬ 
ular clusters (GLC), 2 high mass X-ray binaries (HMB), all 
identified, 295 BL Lac objects (BZBs) 278 flat spectrum ra¬ 
dio quasars (BZQs), 28 non blazar-like active galactic nuclei 
(AGNs) 92 AGN of uncertain type (AGUs) 6 normal galaxies 
and 2 starburst galaxies. The remaining 630 Fermi sources 
are all UGSs. 

A detailed comparison between the IFGL and the 2FGL 
catalog is reported in Nolan et al. (2012). There are 1099 
sources included in both Fermi catalogs with only 381 IFGL 
sources being listed in the IFHL (Ackermann et al. 2013). It 
is worth noting that 267 out of 1099 sources in common be¬ 
tween IFGL and 2FGL changed classification from IFGL 
to 2FGL and among them 172, previously unassociated in 
IFGL, obtained an assigned low-energy counterpart in 2FGL. 
This sample of modified associations includes: 16 sources in¬ 
dicated as generic AGNs that have been classified plus one 
that became a UGS; 54 AGUs that have been recognized 
as blazars (either BZBs or BZQs) and 2 as AGNs; then 17 
blazars actually changed their classification from IFGL to 
2FGL thanks to new spectroscopic data available and, among 
them, 6 actually disappeared since they were listed in the 
2FGL as unassociated as occurred for other 5 sources of 
Galactic origin (i.e., PSRs, PWNe and SNRs). On the other 
hand, there are 172 UGSs of the IFGL that appear to be asso¬ 
ciated in the 2FGL: 20 associated with PSRs, 11 with SNRs 
and 1 with a GLC while all the others with various AGN 
classes (i.e. blazars, radio galaxies, Seyfert galaxies etc.). 

The IFGL catalog includes 216 sources with y-ray anal¬ 
ysis flags and their detection has to be considered care¬ 
fully since they could have been artifacts of the analysis 
(Abdoetal. 2010a). In particular, 103 of these 216 IFGL 
flagged sources have the “c” flag following the IFGL name; 
this indicates that the source is found in a region with 
bright and/or possibly incorrectly modeled diffuse emission. 
On the other hand any non-zero entry in the Flags col¬ 
umn indicates inconsistencies found during the analysis (see 
Abdo et al. 2010a, for details). However we note that 36 out 
of 103 “c”-flagged IFGL sources are also detected in the 
2FGL and only 25 of them were still indicated with a “c” flag 
therein. Thus we decided to keep all the c-flagged sources in 
our analysis since they could be useful for reference in future 
releases of Fermi catalogs. 

3.2. The Fermi-Large Area Telescope Second Source Catalog 

The 2FGL catalog^^ based on 24 months of Fermi operation 
lists 1873 y-ray sources (Nolan et al. 2012). It includes: 

• 108 PSRs, 25 associated plus 83 identified; 

• 3 PWNe, all 3 identified; 

• 61 y-ray sources corresponding to 50 SNRs, 55 associ¬ 
ated plus 6 identified; 

• 11 associated globular clusters (GLC), 4 high mass X- 
ray binaries (HMB), all identified and 1 nova (NOV), 

among the Galactic sources, while in the extragalactic sky: 

• 435 BL Lac objects (BZBs), only 7 identified; 

• 370 flat spectrum radio quasars (BZQs), 17 identified; 

http://fermi.gsfc.nasa.gov/ssc/data/access/lat/2yr_catalog/ 
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• 11 non blazar-like active galactic nuclei (AGNs), 1 
identified; 

• 257 AGNs of uncertain type (AGUs); 

• 6 normal galaxies and 4 starburst galaxies . 

The remaining 577 Fermi sources were are all UGSs. 

The 2FGL catalog lists 472 sources with y-ray analysis 
fiags and 155 of them also have the “c” fiag following the 
2FGL name. It is also worth noting that with respect to the 
IFGL catalog the fraction of UGSs is decreased from 630 out 
of 1451 sources (-43%) to 577 out of 1873 (-31%) in the 
2FGL. 

4. INTRODUCING THE CLASSIFICATION FOR THE LOW-ENERGY 
COUNTERPARTS OF THE Fermi SOURCES 

We adopt the following class definitions for Galactic coun¬ 
terparts of the Fermi sources, summarized in Table 2: high 
mass X-ray binary (hmb), globular cluster (glc), nova (nov), 
millisecond pulsar (msp), pulsar (psr), pulsar wind nebula 
(pwn), supernova remnant (snr), star forming region or Hii 
region (sfr) and binary star (bin). The classes that we as¬ 
sign for extragalactic counterparts are BL Lac object (bzb), 
quasar radio loud of fiat spectrum type (bzq), blazar of un¬ 
certain type (bzu), radio galaxy (rdg), normal galaxy (gal), 
starburst galaxy (sbg) and Seyfert galaxy (sey). We note that 
we used the three-letter designators to indicate the class 
for each counterpart; these are shown in lower case to dif¬ 
ferentiate them from the acronyms used throughout the 
paper (see Table 1). We emphasize that the labels used for 
the extragalactic blazar subclasses are used only if the coun¬ 
terpart corresponds to a source classified in the Roma-BZCAT 
(Massaro et al. 2009; Massaro et al. 2011b). We used ver¬ 
sion 4.1 of the Roma-BZCAT in our analysis. We also em¬ 
phasize that the proposed classification depends on the avail¬ 
ability of the optical spectroscopic information for the low- 
energy counterparts. 

We also considered as a counterpart class blazar candidates 
(ben) to indicate sources having at least one of the following 
requirements: 1) being classified as BL Lac candidate in the 
Roma-BZCAT; 2) sources that belong to the Combined Radio 
All-Sky Targeted Eight-GHz Survey radio catalog (CRATES; 
Healey et al. 2007) with a BL Lac-like or a quasar-like opti¬ 
cal spectrum available in literature; 3) radio sources with a 
blazar-like optical spectrum available in literature (i.e., a pub¬ 
lished spectrum or a description of the features found in the 
near infrared-optical band) and with a WISE counterpart. 

Finally, we labelled as unclassified counterparts (unc) those 
sources satisfying at least one of the following criteria: 1) 
WISE Infrared colors similar to those of y-ray blazars se¬ 
lected according to the association methods developed by 
D’Abrusco et al. (2013) and/or with the kernel density es¬ 
timator analysis (Massaro et al. 2012a; Massaro et al. 2012c, 
see also Section 7 for a refined analysis); 2) a flat radio spec¬ 
trum at frequencies below -1 GHz (i.e., a <0.5); 3) a radio 
source with or without an X-ray counterpart associated with 
the methods adopted in the ILAC and the 2LAC; 4) a radio 
and/or an X-ray source that lies at angular separation between 
the Fermi and the counterpart positions smaller than the max¬ 
imum one for the blazar class (i.e., 0.4525 ° corresponding 
to 2FGL J1746.0-r2316). All the remaining sources are still 
unidentified (ugs). 

The proposed classification scheme has three main advan¬ 
tages. 


1. To have a homogeneous classification of the blazars 
listed in our analysis and to avoid misclassifications we 
adopted the following criterion. We indicate as blazars 
only sources that were previously recognized and clas¬ 
sified in the Roma-BZCAT. There are indeed several 
counterparts classified as BZB or BZQ listed in the 
ILAC and in the 2LAC catalogs that do not have the 
same classification according to the Roma-BZCAT or 
for which no optical spectra were found in literature to 
confirm their classification. 

2. We do not classify a source with a radio and/or X-ray 
counterpart as an active galaxy of uncertain type (AGU) 
since these requirements are not strictly sufficient to 
claim a source as an AGN or as blazar-like and the lack 
of multifrequency information does not allow us to de¬ 
termine precisely its nature. Consequently, we intro¬ 
duced the ben and the unc classes. Once optical spectra 
and multifrequency observations become available a re¬ 
fined classification will be provided. 

3. We introduced the sfr class for cases where Fermi posi¬ 
tion is consistent with a known massive SFR for three 
reasons. First, to highlight the possibility that an un¬ 
known SNR, PWN and/or PSR is embedded therein, 
for which multifrequency observations are necessary 
to confirm their presence. Second, because the use of 
this information could be adopted to refine y-ray dif¬ 
fuse models in future releases of the Fermi catalogs. In 
particular, the largest fraction of these SFR potential as¬ 
sociations were determined using the recent WISE Cat¬ 
alog of Galactic Hii Regions (Anderson et al. 2014), 
which was not available during the preparation of the 
currently-published Fermi catalogs. This could also 
reveal new correspondences between the uncertainties 
of the diffuse y-ray background models and regions of 
interstellar medium in the Galactic plane. Third, this 
could be useful for future investigations of the environ¬ 
ments of the Galactic y-ray sources. 

Finally, we remark that for the sources indicated as candi¬ 
date associations and classified as SFRs, SNRs, PWNs and 
PSRs we investigated the radio and infrared images avail¬ 
able in the NVSS, the Westerbork Synthesis Radio Telescope 
(WSRT) Galactic Plane Compact 327-MHz Source Catalog 
(WSRTGP; Taylor et al. 1996), the Molonglo Galactic Plane 
Survey^^ (MGPS; Green et al. 1999) catalogs and in the WISE 
archives searching for signatures of potential y-ray emitters. 

In Figures 1 and 2 we show two cases of y-ray sources hav¬ 
ing extended structures associated with an open star cluster in 
a SFR found within their Fermi positional uncertainty regions. 
These are two examples of star clusters where some members 
could have been progressed becoming a PSR, PWN, and/or a 
SNR, and so being a plausible counterpart for the y-ray emis¬ 
sion. 

Then in Figures 3 and 4 we report the best two cases of 
SFRs for which a signature of a shell-like or ring-like ex¬ 
tended structures was found within the Fermi positional un¬ 
certainty region that could be ascribed as a shock signature 
potentially responsible for the y-ray emission. 

http://astro.phys.wvu.edu/wise/ 

http://www.astrop.physics.usyd.edu.au/MGPS/ 
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Fig. 1.— The WISE images of the Fermi source lFGLJ0541.1+3542c in 
two different bands at 3.4/im (top), 4.6//m (bottom). In each image the yel¬ 
low ellipse corresponds to the Fermi positional uncertainty region at 68% 
level of confidence. The red contours on the WISE image at 3.4//m indicate 
the radio brightness at 1.4 GHz from the NVSS. The radio-infrared extended 
structure, marked by the green arrow, suggests that a shock occurred in the 
interstellar medium that could be responsible for the y-ray emission. The ra¬ 
dio source is also detected in the Westerbork Northern Sky Survey (WENSS; 
Rengelink et al. 1997) and it is associated with an open star cluster. 

5. CORRELATION WITH EXISTING DATABASES 

We searched the following major radio, infrared, optical 
and X-ray surveys and both the NASA Extragalactic Database 
(NED)^^ and the SIMBAD Astronomical Databaseto verify 
whether multifrequency information can confirm the natures 
of uncertain counterparts and blazar candidates. In the fol¬ 
lowing we list all the catalogs and the surveys searched 
for our investigation; the abbreviations correspond to the 
designators used for the multifrequency notes in Table 4. 

Below ~1 GHz radio frequency we searched the VLA 
Low-Erequency Sky Survey Discrete Source Catalog (VLSS; 
Cohen et al. 2007, - V) and the recent revision VLLSr^^ 

http://ned.ipac.caltech.edu/ 

http://simbad.u-strasbg.fr/simbad/ 

http://heasarc.gsfc.nasa.gov/W3Browse/all/vlssr.html 


Eig. 2.— The WISE images of the Fermi source 1EGLJ1729.1-3641c in two 
different bands at 3.4//m (top), 4.6//m (bottom). In each image the yellow el¬ 
lipse corresponds to the Fermi positional uncertainty region at 68% level of 
confidence. The white radio contours drawn from the NVSS image are over¬ 
laid on the WISE image at 3.4//m. This is an example of a candidate asso¬ 
ciation with a SER that includes the open star cluster 165 from the Northern 
and Equatorial Milky Way catalog built with the Two Micron All-Sky Survey 
(2MASS) observations (Bica et al. 2003). 


(Lane et al. 2014), both the Westerbork Northern Sky Sur¬ 
vey (WENSS; Rengelink et al. 1997, - W) and the Wester¬ 
bork in the Southern Hemisphere Source Catalog (WISH; 
De Breuck et al. 2002, - W), the Sydney University Molon- 
glo Sky Survey (SUMSS; Mauch et al. 2003, - S), the Parkes- 
MIT-NRAO Surveys (PMN; Wright etal. 1994, - Pm), the 
Texas Survey of Radio Sources at 365 MHz (TEXAS; 
Douglas et al. 1996, - T) the Parkes Southern Radio Source 
catalog (PKS; Wright et al. 1990, - Pk), and the Low- 
frequency Radio Catalog of Elat-spectrum Sources (LOR- 
CAT; Massaro et al. 2014a, - L). 

Eor the radio counterparts (i.e., above ~1 GHz), 
we searched the NRAO VLA Sky Survey (NVSS; 
Condon et al. 1998, - N), the VLA Eaint Images of the Ra¬ 
dio Sky at Twenty-Centimeters (EIRST; Becker etal. 1995; 
White et al. 1997, - E), the 87 Green Bank catalog of ra¬ 
dio sources (87GB; Gregory & Condon 1991, - 87) and 
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Fig. 3.— The WISE images of the Fermi source lFGLJ1726.2-3521c in 
two different bands at 4.6//m (top), 12/im (bottom). In each image the yel¬ 
low ellipse corresponds to the Fermi positional uncertainty region at 68% 
level of confidence. The white contours are overlaid on the WISE image at 
llfim are the 1.4 GHz radio emission from the NVSS. The radio-infrared ex¬ 
tended structure, marked by the green arrow, clearly indicates that a shock 
occurred in the interstellar medium within the Fermi positional uncertainty 
region. This could be responsible for the y-ray emission. 

the Green Bank 6-cni Radio Source Catalog (GB6; 
Gregory et al. 1996, - GB), the Australia Telescope 20 
GHz Survey (AT20G; Murphy et al. 2010, - A), the Com¬ 
bined Radio All-Sky Targeted Eight-GHz Survey (CRATES; 
Healey et al. 2007, - c) and the Australian Long Baseline Ar¬ 
ray (LBA) Calibrator Survey of southern compact extragalac- 
tic radio sources^^ (LCSl Petrov et al. 2011, - Ics). 

Eor the infrared, we queried the WISE all-sky survey in the 
AllWISE Source catalog^^ (Wright et al. 2010, - w) and the 
Two Micron All Sky Survey (2MASS; Skrutskie et al. 2006, 
- M) since each WISE source is automatically matched to the 
closest 2MASS potential counterpart (see Cutri et al. 2012, 
for details). 

Then, we also searched for optical counterparts, with or 
http://astrogeo.org/lcs/ 

http://wise2.ipac.caltech.edu/docs/release/allwise/ 


Fig. 4.— The WISE images of the Fermi source 1FGLJ1925.0-I-1720c in 
two different bands at 4.6//m (top), 12//m (bottom). In each image the yellow 
ellipse corresponds to the Fermi positional uncertainty region at 68% level 
of confidence. The white contours are overlaid on the WISE image at 12/im 
are the 1.4 GHz radio emission from the NVSS. The ring-shaped structure, 
indicated by the green arrow, is extended for about 87" and resembles an un¬ 
known SNR or a PWN. It clearly lies within the Fermi positional uncertainty 
region and is also emitting in the radio band. This is a clear example of a 
candidate association with a SFR including a potential SNR/PWN for which 
additional multifrequency observations are necessary to confirm its origin. 

without spectra available, in the Sloan Digitized Sky Sur¬ 
vey Data Release 9 (SDSS DR9; e.g. Ahn et al. 2012, - 
s), in the Six-degree-Eield Galaxy Redshift Survey (6dEGS; 
Jones et al. 2004; Jones et al. 2009, - 6). We also queried 
the United States Naval Observatory (USNO)-Bl Catalog 
(Monet et al. 2003, - U) for the optical counterparts of the ben 
and the unc sources and reported their magnitudes since they 
could be useful to plan future follow up observations. 

In X-rays, we searched the ROSAT all-sky survey in both 
the ROSAT Bright Source Catalog (RBSC; Voges et al. 1999, 
- X) and the ROSAT Eaint Source Catalog (RESC; 
Voges et al. 2000, - X). as well as XMM-Newton Slew Survey 
(XMMSL; Saxton et al. 2008; Warwick et al. 2012, - x), the 
Deep Swift X-Ray Telescope Point Source Catalog (IXSPS; 
Evans et al. 2014, - x) the Chandra Source Catalog (CSC; 
Evans et al. 2010, - x) and Swift X-ray survey for all the Fermi 
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TABLE 2 

Source classification for the y-RAY counterparts. 


TABLE 3 

Multifrequency notes for pulsars. 


I galactic sources 


High mass X-ray binary hmb 

Globural cluster glc 

Nova nov 

Millisecond pulsar (MSP) msp 

Pulsar (PSR) psr 

Pulsar wind nebula (PWN) pwn 

Binary star bin 

Star forming-HII region (SFR) sfr 

(for potential associations with unknown SNRs andjor PWNe and/or PSRs) 
Supernova renmant (SNR) snr 

EXTRAGALACTIC SOURCES 

BE Lac object (BZB) bzb 

Quasar radio loud with flat radio spectrum (BZQ) bzq 

Blazar of uncertain type (BZU) bzu 

Normal galaxy gal 

Radio galaxy rdg 

Starburst galaxy sbg 

Seyfert galaxy sey 

uncertain classification 

Blazar candidate ben 

Unclassified sources unc 

Unidentified gamma-ray source ugs 


UGSs (Stroh et al. 2013). We used the same designator for 
the X-ray catalogs of XMM-Newton, Chandra and Swift. 
These X-ray observatories perform only pointed observa¬ 
tions thus there is the chance that the pointed observa¬ 
tions present in their archives and related to the field of the 
Fermi source was not requested as follow up of the y-ray 
source but for a different reason. Consequently, the dis¬ 
covery of the X-ray counterpart for an associated and/or 
identified source could be serendipitous. 

To perform the searches over all the surveys and the cata¬ 
logs cited above we considered the Icr positional uncertain¬ 
ties reported therein, with the only exceptions being the WISE 
all-sky survey and the SDSS DR9. In these two cases we 
searched the closest IR and optical counterpart within a max¬ 
imum angular separation of 3".3 and of 1".8 for the AllWISE 
survey and the SDSS DR9, respectively. These values have 
been derived on the basis of the statistical analysis described 
in D’Abrusco et al. (2013) and Massaro et al. (2014b) devel¬ 
oped following the approach also presented in Maselli et al. 
(2010) and Stephen et al. (2010). 

To verify the presence of SFRs consistent with the posi¬ 
tional uncertainty of the Fermi sources, we searched in the 
following surveys and catalogs: 1) the Sharpless catalog 
of Hii regions (Sharpless 1959), 2) the Hii regions in the 
(InfraRed Astronomical Satellite (IRAS) point source catalog 
(Hughes & MacLeod; Codella et al. 1994), 3) the radio 
survey of Hii regions at 4.85 GHz (Kuchar & Clark 1997), 
4) the new catalog of Hii regions in the Milky Way 
(Giveon et al. 2005), 5) the radio catalog of Galactic 
HII regions from decimeter to millimeter wavelength 
(Paladini et al. 2003) 6) the SCUBA imaging survey 

of ultra compact Hii regions (Thompson et al. 2006); 
7) the recent WISE Catalog of Galactic Hii Regions 
(Anderson et al. 2014) 8) the QUaD Galactic Plane Survey 
(Culverhouse et al. 2011) 9) the 5 GHz VLA survey of the 
Galacitc plane (Becker etal. 1994) and 10) the WSRTGP 
(Taylor et al. 1996). We remark that for the cross-matching 
performed between the SFR catalogs and the Fermi sources 


Symbol 

Comment 

b 

PSR in a binary system 

e 

PSR detected in y rays by the Compton Gamma-ray Observatory 

g 

PSR discovered in LAT y-ray data 

m 

MSP 

P 

Pulsar discovered by the Pulsar Search Consortium (PSC) 

r 

PSR discovered in the radio 

u 

PSR discovered using a Fermi-hAI seed position 

X 

PSR discovered in X-rays 


we took into account the size of each H ii region in combina¬ 
tion with the Fermi uncertainty position. 

To check for the presence of PSRs we used the Aus¬ 
tralia Telescope National Facility (ATNF) Pulsar Catalog 
(Manchester et al. 2005), the Catalogue of Galactic Super¬ 
nova Remnants (Green 2009), and the census of high-energy 
observations of Galactic SNRs^^ (Ferrand & Safi-Harb 2012). 

Table 4 summarizes all the notes derived from the 
multifrequency analysis. In particular, for the PSRs, 
we report the information shown in the Public List of 
LAT-Detected Gamma-Ray Pulsars^^ and by the Pulsar 
Search Consortium (PSC) as well as from literature (e.g. 

Hobbs et al. 2004; Espinoza et al. 2013; Ray et al. 2012; 

Ray et al. 2013; Pletsch et al. 2012a; Pletsch et al. 2012b; 

Pletsch et al. 2013; Hanabata et al. 2014). Symbols used for 
their multifrequency notes are labelled as reported in Table 3 
(see Abdo et al. 2013, for more details). 

6. THE REFINED ASSOCIATIONS FOR THE Fermi CATALOGS 

The refined associations for the ILGL and the 2LGL catalog 
(hereinafter ILGLR and 2LGLR, respectively) is presented in 
Table 4. Lor each source we report: the ILGL name together 
with the 2LGL and the ILHL ones as derived from the associ¬ 
ations published in those catalogs, the counterpart name, and 
in 37 cases also an alternative association; the category of y- 
ray source association; the class of each counterpart as found 
in the literature; notes derived from the multifrequency inves¬ 
tigation; we also report whether the counterpart lies inside a 
SLR, and specifically for the PSRs in a PWN or a SNR. In ad¬ 
dition, in Table 4 we also report if a PWN or a SNR contains 
a known PSR. The coordinates of each counterpart are also 
included in Table 4 together with the optical magnitudes in B 
and R band from the USNO-Bl Catalog (Monet et al. 2003). 

Within the multifrequency notes in Table 4 we also indicate if 
the spectral energy distribution (SED) of the source is shown 
in Takeuchi et al. (2013), if a redshift estimate is present in 
literature (indicating the reference) and if the radio counter¬ 
part has a fiat radio spectrum (i.e., designated rf in the notes 
of Table 4 whenever a <0.5 in the radio band). In the on-line 
version of the table we also report the infrared analysis fiag 
(Wright et al. 2010) for those sources, classified as blazar-like 
with a counterpart in the AllWISE catalog within 3".3 as well 
as the column with the confirmed redshifts. 

Lor the counterpart names we used the following priorities: 
if the source is a known blazar we adopted the Roma-BZCAT 

http://www.physics.umanitoba.ca/snr/SNRcat/ 

https://confluence.slac.stanford.edu/display/GLAMCOG/Public+List+of+LAT-Detecte 
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nomenclature while for the known pulsars that of the ATNF 
Pulsar Catalog^^ (Manchester et al. 2005). Radio galaxies, 
Seyfert galaxies and starburst names were indicated from 
the Third Cambridge Catalogue of Radio Sources and re¬ 
vised versions (3C and SCR Edge et al. 1959; Bennett 1962; 
Spinrad et al. 1985), the New General catalog and the Messier 
catalog or using their radio names reported in one of the ma¬ 
jor surveys (see Section 5). Then Galactic sources as SNRs 
and PWNe have their names in Galactic coordinates as re¬ 
ported in the Catalogue of Galactic Supernova Remnants^^ 
(Green 2009), while the most common names were used for 
globular clusters and binaries. The counterparts that were 
associated with SFRs were also labelled with their name in 
Galactic coordinates. Finally, the remaining uncertain sources 
were identified by their radio names whenever possible or by 
their ROSAT name one. A handful of counterparts were in¬ 
dicated using the nomenclature of optical catalogs. Adopt¬ 
ing these choices in the names we remark that all the coun¬ 
terparts can be retrieved from the NASA/IPAC Extragalactic 
Database (NED) and SIMBAD archives without using their 
coordinates. 

6.1. The refined associations of the IFGL 

The IFGLR lists 111 identified sources, 880 associated, 
306 candidate associations, and 154 UGSs. The sky distri¬ 
bution of the UGSs is shown in Figure 5 in comparison with 
the IFGL associations in which there were 66 identified ob¬ 
jects, 755 associated, and 630 unidentified. The results of 
our multifrequency analysis for the IFGLR are summarized 
in Table 5 where we list all the sources identified, associ¬ 
ated, and candidate associations for each class separately. We 
remark that the globular cluster GCl 94 (alias M28) associ¬ 
ated with IFGL J1824.5-2449 contains a known MSP PSR 
J1824-2452A while 2MASS-GC01 lies in an SFR. 

Among the PSRs and the MSPs are 21 sources that lie 
within a known PWN and 4 in SNRs as reported in the notes 
of Table 4. Twenty PSRs and MSPs out of 39 indicated as 
candidate associations in both Fermi catalogs have been iden¬ 
tified after the release of the 2FGL catalog, as presented in 
the Public List of FAT-Detected Gamma-Ray Pulsars. This 
strongly supports our introduction of this new category of y- 
ray sources. In five of the 10 associations of PWNe a PSR 
also has been found within the nebula and the same situation 
occurs for 8 SNRs. 

We find that 113 IFGL sources have an SFR within 
their positional uncertainty regions obtained by combining in 
quadrature the y-ray localization at 95% level of confidence 
with the size of the SFR. Three of these sources are the HMBs 
and one is the globular cluster 2MASS-GC01. Fifty-six of 
these 113 have been tentatively associated (i.e., indicated as 
candidate associations) with previously UGSs while the re¬ 
maining Fermi objects all have PSRs, PWNe and/or a SNRs 
embedded in gas clouds and/or interacting with the interstellar 
medium^^. Among the 54 candidate associations with SFRs, 
two contain SNRs, one a PWN, and one a known PSR. How¬ 
ever, the potential associations with SFRs could be used to 
refine the y-ray Galactic diffuse emission models and the as¬ 
sociations in future releases of the Fermi catalogs, as well as 
lead to the discovery of new PWNe and/or unknown SNRs. 

http://www.atnf.csiro.au/research/pulsar/psrcat/ 

http://www.mrao.cam.ac.uk/surveys/snrs/ 

see also http://astronomy.nju.edu.cn/ ygchen/others/bjiang/interSNR6.htm 


None of the MSPs has been found lying within a known SFR 
while this occurs for 21 PSRs, 7 PWNe, and 18 SNRs. 

It is worth noting that IFGL J0503.2-f4526, associated ac¬ 
cording to the 2FGL and the 2LAC analyses with an un¬ 
classified source, has as an alternative potential counterpart a 
SNR that includes a PSR. A similar situation occurs for IFGL 
J1837.5-0659c and IFGL J1846.8-0233c whose y-ray emis¬ 
sion could be ascribed to a PSR in a PWN. In addition, IFGL 
J0622.2-f3751 has as an alternative association with a PSR. 

Within the unclassified candidate associations (i.e., unc), 
we noticed that 8 IFGLR sources are positionally consis¬ 
tent with SFRs. In addition, we also listed in the IFGLR 
the tentative association of IFGL J1653.6-0158 (also known 
as 2FGL J1653.6-0159) with the binary system that includes 
PSR J1653-0158 (Romani et al. 2014). 

Finally, we remark that in the extragalactic sky the 
Small Magellanic Cloud (alias NGC 292) is identified with 
the Fermi source IFGL JOlOl.3-7257 and classified as a 
normal galaxy while the Large Magellanic Cloud is asso¬ 
ciated with all the remaining five Fermi sources listed as 
normal galaxies. 

6.2. The refined associations of the 2FGL 

In the analysis of the IFGLR we have already considered 
the associations for 1099 Fermi sources that are in common 
between the IFGL and the 2FGL, however, we report the 
results for the whole refined 2FGL associations (hereinafter 
2FGLR) in Table 5. As for the IFGLR, the 2FGLR cata¬ 
log is dominated by blazar-like sources in the extragalactic 
sky and by pulsars around the Galactic plane. The 2FGLR 
lists 126 identified sources as in the original catalog, 1176 
associated, 282 candidate associations and 289 UGSs. The 
2FGLR sources are classified as: 4 high mass X-ray bi¬ 
naries, 11 globular clusters, 51 MSPs, and 98 PSRs, 21 that 
lie in a PWN and 6 in an SNR, and 17 associated with a 
SFR indicated in our Table 4. In the case of 2FGLR, among 
the PSR candidate associations that are not in the IFGLR 
only one source has been recently identified according to the 
Public List of FAT-Detected Gamma-Ray Pulsars (see also 
Abdo et al. 2013). There are also 9 PWNe and 66 sources 
associated with SNRs, 11 of them with a known PSR lying 
inside the remnant. Then we list 63 Fermi sources as candi¬ 
date associations with SFRs, one also including a SNR. With 
respect to the IFGL there is also a new association with a nova 
(e.g., Cheung et al. 2014). 

The extragalactic sky includes 277 BZBs, 338 BZQs, 54 
BZUs, 3 Seyfert galaxies, and 3 starburst galaxies. Four Fermi 
sources are associated with the y-ray emission arising from 
the Large Magellanic cloud, as occurred for the objects clas¬ 
sified as galaxy in the IFGLR, and one with M31 (i.e., the 
Andromeda Galaxy) (Abdo et al. 2010d). In addition we find 
271 blazar candidates (ben), 320 Fermi unclassified sources 
and 289 UGSs. In Figure 6 we report the comparison be¬ 
tween the sky distributions of all the UGSs previously listed 
in the 2FGL and the remaining ones determined in this refined 
association analysis. 

7. INFRARED COLORS OF y-RAY BLAZAR CANDIDATES AND 
UNCLASSIFIED SOURCES 

As already done for the AGUs (Massaro et al. 2012a) and 
the UGSs (e.g., Paggi et al. 2013) listed in the 2FGL we per¬ 
formed a non-parametric analysis of the infrared colors for the 
Fermi sources classified as ben and unc in the merged Fermi 
catalog IFGLR-I-2FGLR. We used kernel density estimation 
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TABLE 4 

Refined associations of the Fermi catalogs (first 30 lines). 



Column description. (1): IFGL name; (2): 2FGL name; (3): IFHL name; (4): y-ray source category (I=identified; A=associated; C=candidate; U=unidentified; 
see Section 2 for details); (5): name of the first low-energy counterpart; (6): multifrequency notes on the first low-energy counterpart; (7): redshift estimate if the 
source is extragalactic (in the on line version of the catalog there is also a note if the z measurement is uncertain); (8): notes if the first low-energy counterpart 
lies inside a SFR/SNR/PWN; (9): notes if the first low-energy counterpart includes a PSR/MSP; (10): source class for the first low-energy counterpart (see 
Section 4 for details); (11): R.A. (J2000) for the first low-energy counterpart; (12): Dec. (J2000) for the first low-energy counterpart; (13): B magnitude for 
the first low-energy counterpart taken from the USNO-Bl Catalog (Monet et al. 2003), only for the ben and the unc sources; (14): R magnitude for the first 
low-energy counterpart taken from the USNO-Bl Catalog (Monet et al. 2003), only for the ben and the unc sources; (15): name of the second low-energy 
counterpart; (16): multifrequency notes on the second low-energy counterpart; (17): source class for the second low-energy counterpart (see Section 4 for 
details); (18): R.A. (J2000) for the second low-energy counterpart; (19): Dec. (J2000) for the second low-energy counterpart; (20): B magnitude for the 
second low-energy counterpart taken from the USNO-Bl Catalog (Monet et al. 2003), only for the ben and the unc sources; (21): R magnitude for the second 
low-energy counterpart taken from the USNO-Bl Catalog (Monet et al. 2003), only for the ben and the unc sources; (22): probability that the first low-energy 
counterpart has infrared colors similar to the known IFGL BZBs estimated with the KDE method (see Section 7 for details); (23): probability that the first 
low-energy counterpart has infrared colors similar to the known IFGL BZQs estimated with the KDE method (see Section 7 for details). 

Symbols used for the multifrequency notes are all reported in Section 5 together with the references of the catalogs/surveys. 

References for the optical spectra reported in the table: Ackermann et al. (2011), Ahn et al. (2012), Baker et al. (1995), Bauer et al. (2000), Bade et al. (1995), 
Beuermann et al. (1999), Bikmaev et al. (2008), Britzen et al. (2007), Drinkwater et al. (1997), Healey et al. (2008), Hewett & Wild (2010), Hewitt & Burbidge 
(1989), Johnston et al. (1995), Jones et al. (2004), Jones et al. (2009), Landoni et al. (2014) Landt et al. (2001), Lister et al. (2011), Mao (2011), Maza et al. 
(1995), Marti et al. (2004), Masetti et al. (2013), Massaro et al. (2014c) Mitton et al. (1977), Paggi et al. (2013), Quintana & Ramirez (1995), Schwope et al. 
(2000), Shaw et al. (2013a), Shaw et al. (2013b), Stern & Assef (2013), Shekel & Kuehr (1996a,b) Titov et al. (2011), Vandenbroucke et al. (2010), Vettolani et 
al. (1989), Takeuchi et al. (2013), Thompson et al. (1990), Tsarevsky ert al. (2005), White et al. (1988). 


IFGL unidentified gamma-ray sources 


IFGLR unidentified gamma-ray sources 




Fig. 5. — The all-sky distribution of the UGSs in the IFGL (left) in comparison with those listed in the IFGLR (right panel). Red circles mark those UGSs 
with no y-ray analysis flags and blue squares those with a flag. The projection is Hammer-Aitoff projection in Galactic coordinates 


(KDE; see e.g., D’Abmsco et al. 2009; Laurinoetal. 2011, 
and reference therein) to verify the consistency of their 
infrared colors with the so-called WISE Gamma-ray 
Strip (e.g., Massaro et al. 2011a; Massaro et al. 2012a; 
Massaro et al. 2013c) The KDE technique allows us to 
estimate the probability function of a multivariate distribution 
and does not require any assumption about the shape of the 
“parent” distributions. 

Consequently, for all the ben and unc sources with a WISE 
counterpart we can provide an estimate of the probability 
Tikde that such source is consistent with the WISE Gamma- 
ray Strip. In particular we differentiated between nue,hzh 


and TTkdeMq considering the comparison with the infrared 
colors of the BZB and BZQ subclasses, respectively. We 
considered as training samples to build the WISE Gamma- 
ray Strip(Massaro et al. 2011a; Massaro et al. 2013a) and to 
compare the infrared colors, the Roma-BZCAT blazars listed 
in both IFGLR and 2FGLR. We remark that WISE counter¬ 
parts of the ben and of the unc sources with the infrared anal¬ 
ysis flags^^ (Wright et al. 2010) have been considered in our 

As for example contamination and con¬ 
fusion from nearby bright sources, see e.g., 
http://wise2.ipac.caltech.edu/docs/release/allsky/expsup/sec6_3a.html 
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2FGL unidentified gamma-ray sources 


2FGLR unidentified gamma-ray sources 




Fig. 6.— All-sky distribution of the UGSs in the 2FGL (left) in comparison with those listed in the 2FGLR (right panel). Red circles mark those UGSs with no 
y-ray analysis flags and blue squares those with a flag. 


TABLE 5 

Summary of categories eor the y-RAY source associations in the reeined 
ASSOCIATION LIST OE BOTH Fermi CATALOGS. 


1 1 

1 

lEGLR 


1 

1 

2EGLR 


1 

1 Class 

1 ID- 

ASS. 

CAN. 

Tot. 

1 ID- 

ASS. 

CAN. 

Tot. 1 

GALACTIC 1 1 

bin 

0 

1 

0 

1 

0 

1 

0 

1 

hmb 

3 

0 

1 

4 

4 

0 

0 

4 

glc 

0 

9 

0 

9 

0 

11 

0 

11 

nov 

0 

0 

0 

0 

1 

0 

0 

1 

msp 

26 

16 

8 

50 

26 

17 

8 

51 

psr 

53 

8 

31 

92 

56 

12 

30 

98 

pwn 

0 

8 

2 

10 

2 

6 

1 

9 

snr 

3 

38 

5 

46 

6 

55 

5 

66 

sfr 

0 

0 

54 

54 

0 

0 

63 

63 

EXTRAGALACTIC 



bzb 

7 

231 

1 

239 

7 

270 

0 

277 

bzq 

15 

271 

1 

287 

18 

320 

0 

338 

bzu 

2 

46 

0 

48 

2 

52 

0 

54 

gal 

1 

5 

0 

6 

2 

4 

0 

6 

rdg 

0 

5 

0 

5 

1 

7 

0 

8 

sbg 

0 

3 

0 

3 

0 

3 

0 

3 

sey 

0 

2 

0 

2 

0 

3 

0 

3 

UNCERTAIN 



ben 

0 

188 

25 

213 

0 

242 

29 

271 

unc 

1 

49 

178 

228 

1 

173 

146 

320 

ugs 

0 

0 

0 

154 

0 

0 

0 

289 


Col. (1) Source class. Col. (2) ID.= identified. Col. (3) ASS.= associated. Col. (4) CAN.= candidate associations. 

analysis and their fiags are reported in Table 4 together with 
the probabilities derived from the KDE analysis. 

In Figure 7, the isodensity contours drawn from the KDE 
density probabilities are plotted for the E.^rm/-Roma-BZCAT 
blazars in the [3.4]-[4.6]-[12] yum color-color diagram. The 
IR colors of the sources listed in the refined association list of 
the Fermi catalogs classified as ben and unc are also shown. 
Values of TikdeMb ^kde,bzq than 5% indicate that the 

source has infrared colors consistent with the portion WISE 
Gamma-ray Strip constituted by the BZBs and by the BZQs, 
respectively, at 95% confidence level. 

Our analysis confirms that 419 out of 499 total ben (262) 
and unc (237) sources in the refined and merged associa¬ 
tions list of the Fermi catalogs are consistent with the WISE 

for additional details 


Gamma-ray Strip at 95% confidence level. For the unc 
sources this situation occurs for 172 sources out of 237 ob¬ 
jects while for the ben in 242 out of 262 cases. We remark 
that 374 out of the 499 sources analyzed with the KDE are as¬ 
sociated, 1 is identified, and the remaining 124 are candidate 
associations. This analysis supports the introduction of the 
category for the candidate associations, since a large fraction 
of the ben and the unc with a WISE counterpart are likely to be 
blazar-like sources being consistent with the WISE Gamma- 
ray Strip. 

In our analysis, we took into account the correction for 
Galactic extinction for all the WISE magnitudes according 
to the Draine (2003) relation. As shown in D’Abrusco et al. 
(2013), this correction affects only marginally the [3.4]-[4.6] 
yum color, mostly at low Galactic latitudes (i.e., \b\ <15 °). 

8. OPTICAL SPECTROSCOPIC OBSERVATIONS OF y-RAY BLAZAR 
CANDIDATES 

The total number of Fermi sources listed uniquely in the 
IFGLR and in the 2FGLR are 2219. We found spectro¬ 
scopic information for 177 of them not reported in the pre¬ 
vious IFGL and 2FGL catalogs, including 8 observed during 
our spectroscopic campaign (see below). 

We have been able to confirm spectroscopically 117 BE 
Lacs, 21 with firm redshift estimates and 27 quasars (QSOs); 
an additional eight sources appear to have a BE Lac nature but 
the lack of good optical spectra in the literature or, if a spec¬ 
trum is available, a low signal to noise ratio did not allow us 
to verify their natures. Thus 24 out of these 177 sources still 
have an uncertain classification and an uncertain redshift esti¬ 
mate. In particular among them we found 11 sources listed as 
blazars in Ackermann et al. (2011) plus one listed in Healey 
et al. (2007) with no optical spectra published or described in 
literature. These 177 sources are all listed in Table 4 as ben’ 
with the only exception being one radio source that appears 
to be a radio galaxy; in the table notes we report their possi¬ 
ble classifications and their redshift estimates with a question 
mark (?) indicating that further investigation is required. 

Optical spectroscopic observations of four candidate as¬ 
sociations and of three associated with unclassified sources 
(unc) were performed with the 2.1-m telescope of the Obser- 
vatorio Astronomico Nacional (OAN) in San Pedro Martir 
(Mexico) on the nights between 28 June and 2 July 2014 
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Fig. 7.— The isodensity contours drawn from the KDE technique for the 
F<?rmi-Roma-BZCAT blazars (i.e., training sample) that constitute the WISE 
Gamma-ray Strip (Massaro et al. 2011a) in the [3.4]-[4.6]-[12] jim color- 
color plot. Blue contours were computed from the infrared colors of the 
BZBs while the red ones were derived from the BZQs. The infrared colors 
of the sources classified as ben (top - black circles) and unc (bottom - green 
squares) are also reported to show their consistency with the WISE Gamma- 
ray Strip. The numbers appearing close to each contour correspond to the 
values of nkde in both panels. The last contour correspond to 95% confidence 
level (see Section 7 for details). A similar analysis for the active galaxies of 
uncertain type listed in the 2FGL was aheady performed in Massaro et al. 
(2012a). Here this study has been repeated for both the IFGLR and 2FGLR 
lists, refined association list of the Fermi catalogs. 

(UT) as reported in Table 6. The telescope carries a Boiler 
& Chivens spectrograph and a 1024x1024 pixel E2V-4240 
CCD. The spectrograph was tuned in the 4000-^8000 A range 
(grating 300 l/mm), with a resolution of 4.5 A per pixel, which 
corresponds to 8 A (Full-width-half-maximum), and a 2'.'5 
slit. Data were wavelength calibrated using Copper-Helium- 
Neon-Argon lamps, while for flux calibration spectrophoto- 
metric standard stars were observed twice during every night 
of the observing run. 

In addition we also observed the y-ray blazar candidate 
BZB J2340-I-8015 associated with IFGL J2341.6-1-8015 at the 
Observatorio Astroflsico Guillermo Haro (OAGH) located at 
Cananea, Sonora in Mexico. The telescope was equipped with 
a spectrograph Boiler & Chivens with a CCD SITe Ikxlk, 


tuned in the 4000-^7000 A range (grating 150 l/mm). We used 
a slit width of 275 with a resolution ~15A (FWHM). 

The data reduction for both telescopes and instruments was 
carried out using the Image Reduction and Analysis Facility 
(IRAF) package developed by the National Optical Astron¬ 
omy Observatory, including bias subtraction, spectroscopic 
flat Adding, optimal extraction of the spectra and interpola¬ 
tion of the wavelength solution. All spectra were reduced 
and calibrated employing standard techniques in IRAF and 
our own IDL routines (see also Matheson et al. 2008, for ad¬ 
ditional details). 

These eight sources for which new spectroscopic data are 
provided in our analysis are listed in Table 4 as y-ray blazar 
candidates (i.e., ben). Five of them belong to the IFGLR and 
are all classifled as BL Lac objects according to our analy¬ 
sis; for three of them the presence of several spectral features 
allowed us to determine their redshifts. The log of the spectro¬ 
scopic observations and the results of our analysis are reported 
in Table 6, while in Figure 8, Figure 9, Figure 10 and Fig¬ 
ure 11 we show their optical spectra together with their flnd- 
ing charts. The remaining three sources are two quasars cor¬ 
responding to 2FGL J1848.6-f3241 and 2FGL J202L5-f0632 
as shown in Figures 12 and 13 at redshifts 0.981 and 0.217, 
respectively, and 2FGL J203L0-f 1938, a BL Lac object with 
an uncertain redshift since there is a unique emission line vis¬ 
ible, potentially identifled as Mg ii (see Figure 14). 

During our observing nights we also obtained the spec¬ 
tra for two 2FGL sources, 2FGL J1719.3-1-1744 and 2FGL 
J180L7-f4405, each with uncertain redshift as reported in 
the Roma-BZCAT (see the flgures in the Appendix). We 
were able to confirm the redshift for 2FGL J1801.7-1-4405 
while 2FGL J1719.3-f 1744 (alias BZB J1719-f 1745) is com¬ 
pletely featureless. Results for these two additional spec¬ 
tra are also reported in Table 6. We also observed IFGL 
J1942.7-f 1033, which is a BL Lac object already classifled in 
literature (Tsarevsky et al. 2005; Masetti et al. 2013, and Ap¬ 
pendix for the figure), and IFGL J2300.4-r3138 for which a 
series of unidentified absorption features are clearly visible 
in its optical spectrum. Assuming that these unidentified ab¬ 
sorption features are due to Mg ii intervening systems along 
the line of sight, the source should lie at redshift >0.96. This 
source was in the sample observed by Shaw et al. (2013a) 
but we cannot confirm their results since we did not And any 
spectral feature identifiable as Civ (see Figure 15). Finally 
we also report the spectrum of BZB J2340-r8015 associated 
with IFGL J234L6-I-8015 observed at OAGH that was orig¬ 
inally classifled as a BL Lac candidate in the Roma-BZCAT 
at redshift 0.274. We have been able to confirm the BL Lac 
nature of this source but not its redshift due to its featureless 
spectrum (see Figure 16). All our finding charts are from the 
Digitized Sky Survey^^. 

9. COMPARISON WITH STATISTICAL ANALYSES 

We compared our multifrequency analysis on the candidate 
associations with the results obtained via statistical analyses 
based on the Classification Tree (CT) and on the Logistic 
Regression (LR) procedures performed on the IFGL sources 
(Ackermann et al. 2012) and using the Random Forest algo¬ 
rithm, named sybil, presented by Mirabal et al. (2012) for 
the 2FGL sources. 

http://archive.eso.org/dss/dss 

http://www.gae.ucm.es/ mirabal/sibyl.html 
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Fig. 8.— Upper panel) The optical spectra of the counterpart associated with 
IFGL J1548.7+6311 observed at OAN in San Pedro Martir (Mexico) on 28 
June 2014. The absorption features identified as Can H&K and the G band 
used to determine its redshift are marked. The source has been classified as 
a BL Lac. (Lower panel) The 5' x 5' finding chart from the Digitized Sky 
Survey (red filter). The potential counterpart of IFGL J1548.7+6311, the 
target of our observation, is indicated by the red circle. 


TABLE 6 

Log and results of optical spectroscopic observations. 


Fermi 

Obs. Date 

Exposure 

Class 

z 

Spectral 

Name 

yyyy-mm-dd 

(sec) 



Features 

1FGLJ1548.7+6311 

2014-06-28 

2x1800 

BL Lac 

0.269 

Can H&K, G, H;S, Mgl 

1FGLJ1844.1+1547 

2014-06-29 

2x1800 

BL Lac 

? 

none 

IFGL J2014.4+0647 

2014-06-29 

2x1800 

BL Lac 

0.341 

Can H&K, G, Hj8, Mgl 

1FGLJ2133.4+2532 

2014-06-29 

2x1800 

BL Lac 

0.294 

Can H&K, G, Hj8, Mgl 

2FGLJ1848.6+3241 

2014-07-02 

2x1800 

QSO 

0.981 

Mgn, [Oil] 43727 

2FGLJ2021.5+0632 

2014-07-02 

1800 

QSO 

0.217 

Hy, HyS, [OIII] 44959 & 45007 

2FGLJ2031.0+1938 

2014-07-01 

2x1800 

BL Lac 

0.668? 

Mgn? 

2FGLJ1719.3+1744 

2014-06-30 

2x1800 

BL Lac 

? 

none 

2FGLJ1801.7+4405 

2014-06-30 

2x1800 

QSO 

0.663 

Mg II, [Oil] T3727 

1FGLJ1942.7+1033 

2014-06-29 

2x1800 

BL Lac 

? 

none 

IFGL J2300.4+3138 

2014-07-01 

2x1800 

BL Lac 

? 

unid. abs. systems 

1FGLJ234L6+8015 

2014-08-28 

3x1800 

BL Lac 

? 

none 


For the IFGL sources the analysis executed with both the 
CT and the LR statistical methods assesses the probability 
of correct classification based on fitting a model form to the 
Fermi data. The result is a PSR-like or an AGN-like classifi¬ 
cation, occurring when the y-ray source properties are likely 
more consistent with those of a pulsar or an active galaxy, 
respectively; a similar prediction is also made for the 2FGL 
sources analyzed by the sybil algorithm. Thus, all these meth- 


Fig. 9.— Upper panel) The optical spectra of the counterpart associated 
with IFGL J1844.1+1547 observed at OAN in San Pedro Martir (Mexico) 
on 29 June 2014. The source has been classified as a BL Lac on the ba¬ 
sis of its featureless continuum. (Lower panel) The 5' x 5' finding chart 
from the Digitized Sky Survey (red filter). The potential counterpart of IFGL 
J1844.1 + 1547, the target of our observation, is indicated by the red circle. 

ods provide a classification of the potential counterpart on the 
basis of the y-ray properties (e.g., spectral shape, variability, 
etc.) but cannot permit to locate it. Consequently we can only 
verify if a candidate association that appear to be a blazar or a 
pulsar have the corresponding y-ray source classified as AGN- 
like or PSR-like by the statistical analyses cited above. 

There are 286 candidate associations of type bzb in the 
IFGLR that were analyzed in Ackermann et al. (2012). 
Within this sample we found that there is only one bzb and 
one hmb, the first statistically classified as AGN-like while 
the second as a PSR-like. Then there are 24 sources classi¬ 
fied as ben, with 23 of them indicated as AGN-like while 
the last is classified as PSR-like. 

Among the Galactic sources, there are: 

• 7 sources, classified by us as potential MSPs, four ap¬ 
pear to be PSR-like while 3 show y-ray properties more 
similar to AGNs (i.e., AGN-like). 

• 25 out of 29 Galactic sources classified are in our analy¬ 
sis as PSRs are confirmed PSR-like by statistical meth¬ 
ods while 4 are indicated as AGN-like; 

• 2 PWN candidate associations, both resembling PSR- 
like sources in y-rays. 

• 5 sources reported in our Table 4 as SNRs, four are clas¬ 
sified as PSR-like objects and one has an uncertain sta- 
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IFGL J2014.4+0647 


IFGL J2133.4+2532 
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Fig. 10.— upper panel) The optical spectra of the counterpart associated 
with IFGL J2014.4+0647 observed at OAN in San Pedro Martir (Mexico). 
The source was observed twice on 29 June 2014 and on 30 June 2014 and both 
spectra are shown. The absorption features identified as Ca ii H&K and the G 
band used to determine its redshift are marked. The source has been classified 
as a BL Lac. Lower panel) The 5' x 5' finding chart from the Digitized 
Sky Survey (red filter). The potential counterpart of IFGL J2014.4+0647 is 
indicated by the red circle. 

tistical classification (i.e., indicated as conflict in Ack- 
ermann et al. 2012). 

• 33 out of 47 candidate associations classified as SFRs 
are PSR-like while only 9 appear to be classifiable as 
AGN-like by the statistical methods; the remaining five 
are unclassified. 

The comparison between our candidate associations in the 
Galactic plane and the results of the statistical analyses sup¬ 
ports the relevance of searching for the presence of SFRs 
within the Fermi positional uncertainty region where PSRs, 
SNRs, PWNe could be embedded therein or shocked regions 
could be potential sources of y-rays not yet confirmed. Fi¬ 
nally, in the sample of 170 unc sources, classified as candi¬ 
date associations by our analysis, 125 of them appearing to be 
AGN-like while 40 are indeed PSR-like; the remaining five 
are conflicts between the LR and the CT methods. In partic¬ 
ular, four out of eight unc sources having a SFR consistent 
with the Fermi position are classified as PSR-like on the basis 
of their y-ray behavior. 

In the case of the 2FGLR, we investigated 133 candidate 
associations with the sybil procedure, and classified them ac¬ 
cordingly. Of these, in the sample of 107 Fermi sources that 
are classified as AGN-like by sybil. 24 are ben, 78 are unc 


4U.U-# ^ j2133.4+25>2 

* 

31:00.0 ' A 


2''9:00.0 * • • (*) • 

• * 

^ 0.0 

28:00.0 * ' ■* ' •• 



27:00.0 . * • 

22.0 _ 16.0- 21:*3:10.0 04.0 

Fig. 11. — upper panel) The optical spectra of the counterpart associated 
with IFGL J2133.4+2532 observed at OAN in San Pedro Martir (Mexico) 
on 29 June 2014. The absorption features identified as Can H&K and the G 
band used to determine its redshift are marked. The source has been classified 
as a BL Lac. Lower panel) The 5' x 5' finding chart from the Digitized 
Sky Survey (red filter). The potential counterpart of IFGL J2133.4+2532 is 
indicated by the red circle. 

type; then the remaining objects are two SFRs, three PSRs, 
seven MSPs (three expected to be PSR-like y-ray sources) and 
five SFRs all classified as AGN-like. 

Finally, we highlight that the results of the statistical anal¬ 
yses are in agreement with the classification proposed in the 
refined association lists of both Fermi catalogs. 

10. GAMMA-RAY CONNECTIONS 
10.1. The radio-y-ray connection 

Many attempts were made in the past to test 
for correlations between the radio and y-ray emis¬ 
sions of AGNs and in particular for blazars 
(e.g.. Sleeker etal. 1993; Padovani et al. 1993; 

Salamon & Stecker 1994; Taylor et al. 2007). This connec¬ 
tion was also used before the Fermi era to search for counter¬ 
parts of the y-ray sources (Mattox et al. 1997), so motivating 
all the past and present radio follow up campaigns for the 
Fermi sources (e.g., Petrov et al. 2013; Schinzel et al. 2014). 

However, biases and selection effects have to be taken into 
account to prove this correlation properly, since it is important 
to address intrinsic source variability, biases due to redshift 
dependence (Elvis et al. 1978), the “common distance” bias 
(Pavlidou et al 2011), problems related to source misidentifi- 
cations and incorrect associations, to name a few. All these 
issues can mimic a correlation (Mucke et al. 1997). 
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Fig. 12.— Upper panel) The optical spectra of the counterparts associated 
with 2FGL J1848.6+3241 observed at OAN in San Pedro Martir (Mexico) on 
2 July 2014. The emission lines marked allowed us to estimate their redshift. 
The source has been classified as a QSO. Lower panel) The 5' x 5' finding 
chart from the Digitized Sky Survey (red filter). The potential counterpart of 
2FGL J1848.6+3241 is indicated by the red circle. 

Since the launch of Fermi several investigations have 
also been performed to search for a definitive an¬ 
swer regarding the existence of the radio-y-ray con¬ 
nection (Ghirlanda et al. 2010; Ghirlanda et al. 2011; 
Mahony et al. 2010) until it was proved and described 
accurately in Ackermann et al. (2011b). Recently we 
also suggested that a link between the radio and the 
y-ray emissions in blazars can be extended well be¬ 
low ~1 GHz (Massaro et al. 2013b; Nori et al. 2014; 
Massaro et al. 2013d). 

Given the new multifrequency analysis carried out here 
and the candidate associations, in particular for the ben class 
we illustrate the current status of the radio-y-ray connection. 
To this end we show both the fiux-fiux and the luminosity- 
luminosity scatter plots for the blazars in our merged refined 
list of Fermi associations (see Figure 17 and Figure 18). 

A study of the correlations is out of the scope of this pa¬ 
per. Our main goal is verify that sources classified as ben 
are in agreement with the behavior of blazars in these param¬ 
eter spaces. It is worth noting that we divided our sample 
of blazar-like sources in two subsamples: the northern one 
and the southern one. the former includes all the blazars (i.e., 
bzb, bzq, and ben classes) that have a radio counterpart within 
the NVSS footprint while in the latter the subsample includes 
those in the area covered by the SUMSS. We distinguish these 
two samples since the radio surveys were carried out at differ- 


Fig. 13. — Upper panel) The optical spectra of the counterparts associated 
with 2FGL J2021.5+0632 observed at OAN in San Pedro Martir (Mexico) on 
2 July 2014. The emission lines marked allowed us to estimate their redshift. 
The source has been classified as a QSO. Lower panel) The 5' x 5' finding 
chart from the Digitized Sky Survey (red filter). The potential counterpart of 
2FGL J2021.5+0632 is indicated by the red circle. 

ent frequencies. 

Finally, we remark that only blazars with a firm redshift es¬ 
timate as reported in our merged list of refined associations 
are shown in Figure 18. According to this figure we could ex¬ 
pect most of the ben sources to be BL Lac objects, although 
this has to be confirmed with optical spectroscopic observa¬ 
tions. 

10.2. The infrared-y-ray connection 

D’Abrusco et al. (2012) studied the Fermi blazar sample 
listed in the 2FGL and lying in the footprint of the WISE Pre¬ 
liminary survey and reported the discovery of a correlation 
between their infrared and their y-ray spectral indices. This is 
directly related to their peculiar infrared colors and the WISE 
Gamma-ray Strip (Massaro et al. 2011a). 

Here we report an updated scatterplot for the spectral in¬ 
dices (see Figure 19) as originally presented in D’Abrusco 
et al. (2012). The correlation the spectral shapes in the in¬ 
frared and in y-rays is still present and as occurred in the case 
of the radio - y-ray scatter plot (see Section 10.1) the loca¬ 
tions of the ben sources appear to be more consistent with 
the BL Lac population. The linear correlation coefficient for 
the whole data set is 0.64. Finally, we also present the con¬ 
nection between infrared and y-ray fluxes and luminosities in 
Figure 20 that is again in good agreement with our previous 
results (D’Abrusco et al. 2012; D’Abrusco et al. 2013). 
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Fig. 14.— Upper panel) The optical spectra of the counterpart associated 
with 2FGL J2031.0+1938 observed at OAN in San Pedro Martir (Mexico) 
on 1 July 2014. The emission line potentially identified as Mgi is marked. 
The source has been classified as a BL Lac. Lower panel) The 5' x 5' finding 
chart from the Digitized Sky Survey (red filter). The potential counterpart of 
2FGL J2031.0+1938 is indicated by the red circle. 

11. SUMMARY AND CONCLUSIONS 

Two years after the release of the second Fermi source cata¬ 
log we present a comprehensive multifrequency investigation 
of all the y-ray associations listed in both the IFGL and the 
2FGL catalogs. 

First, we introduced a new category of y-ray source associ¬ 
ations in addition to identified and associated sources. We 
label as candidate associations those y-ray sources having 
a potential low-energy counterpart of few specific classes of 
well-known y-ray emitters lying within the Fermi positional 
uncertainty region and/or with angular separations between 
the Fermi and the counterpart position smaller than the max¬ 
imum one for all the associated sources of the same class. 
Then we provided a new classification scheme for the low- 
energy counterparts of the Fermi sources based on the multi¬ 
frequency observations and on the optical spectroscopic in¬ 
formation now available. We also presented here a cross¬ 
matching between the Fermi catalogs and several surveys of 
SFRs (see Section 5 for details) to highlight the possibility 
that an unknown SNR, PWN and/or PSR is embedded therein, 
and to provide information that could be used to refine mod¬ 
els of diffuse Galactic y-ray for future releases of the Fermi 
catalogs. 

The total number of y-ray sources considered in our in¬ 
vestigation for both the IFGL and the 2FGL comprises 2219 
unique Fermi objects, all listed in Table 4 with each assigned 


Fig. 15.— Upper panel) The optical spectra of the counterparts associated 
with 2FGL J2300.4+3138 observed at OAN in San Pedro Martir (Mexico) 
on 1 July 2014. Several unidentified absorption features superimposed to 
the optical continuum are visible; however, we were not able to identify the 
Civfeature previously reported by Shaw et al. (2013a). (Lower panel) The 
5' X 5' finding chart from the Digitized Sky Survey (red filter). The poten¬ 
tial counterpart of IFGL J2300.4+3138 pointed during our observations is 
indicated by the red circle. 

counterpart and their main multifrequency properties. Over¬ 
all, in the refined association list of the Fermi catalogs, we 
found 174 Fermi sources with a SFR consistent with their y- 
ray positions. In particular, 60 Fermi objects out of these 174 
do not have the “c” flag in at least one of the Fermi names 
and include: 13 identified sources, 17 associated, and 30 can¬ 
didate associations. Their counterparts are also classified as 
four HMBs; 17 PSRs, among which three lie in a PWN and 
two in a SNR; three PWNe, one with a PSR; one radio galaxy; 
thirteen SNRs, two with a PSR included; one binary star; three 
unclassified sources; and 19 SFRs. 

We found spectroscopic information for 177 extragalactic 
y-ray sources not reported in the previous version of the Fermi 
catalogs. We included analyses of eight new optical spectro¬ 
scopic observations performed with the 2.1-meter telescope of 
the OAN in San Pedro Martir and with the OAGH to confirm 
the natures of these blazar-like sources found in our multi¬ 
frequency investigation. Blazar candidates and counterparts 
with unknown origin (i.e., unc) were also analyzed with the 
KDE technique to determine the fraction having WISE IR col¬ 
ors similar to known y-ray blazars. We compared the refined 
associations listed in both the Fermi catalogs with the clas¬ 
sification proposed on the basis of statistical investigations 
(Ackermann et al. 2012; Mirabal et al. 2012) and we found a 
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IFGL J2341.6+8015 


Radio - y-ray flux densities: Northern sample 



2FGL J2341.6+8015 


30.0 23:40:00 


Fig. 16.— Upper panel) The optical spectra of the counterparts associated 
with IFGL J2341.6+8015 observed at OAGH on 28 August 2014. The source 
has been classified as a BL Lac on the basis of its featureless continuum. 
Lower panel) The 5' x 5' finding chart from the Digitized Sky Survey (red 
filter). The potential counterpart of IFGL J2341.6+8015 is indicated by the 
red circle. 


good agreement with our results. 

We note that in the combined list of refined associations for 
IFGL and 2FGL, among 2219 sources, 394 are still UGSs 
(i.e., ~ 18% of the entire sample) with 191 of them having no 
y-ray analysis fiags. This clean sample of 191 UGSs appears 
to have a uniform distribution in the sky with a small excess 
toward the Galactic plane, as shown in Figure 21. Moreover, 
we conclude that the fraction of Fermi sources with plausible 
counterparts, combining identifications, associations and can¬ 
didate associations, within their y-ray positional uncertainty 
regions, is -80% and up to -90% when considering sources 
with no y-ray analysis fiags. 

Finally, we remark that within the sample of 191 UGSs that 
do not have any y-ray analysis fiags are 40 sources (listed in 
Table 11) that do not have any NVSS and/or SUMSS radio 
source within their positional uncertainty regions at 95% level 
of confidence. Their all-sky distribution, separated into those 
lying in the NVSS footprint (i.e., declination greater than -40 
°) and the others observable from the Southern hemisphere is 
shown in Figure 22. The excess in the subsample of those 
with declinations less than -40°, clearly visible in Figure 22, 
could be due to shallower coverage by SUMSS relative to the 
NVSS catalog. It is worth mentioning that for these 40 UGSs 
the distributions of the y-ray spectral index and of the y-ray 
fiuxes do not show any significant differences from those of 
BZQs and the BZBs (see Figure 23). However it is intriguing 


Radio - Y-ray flux densities: Southern sample 



So.843GHz("*Jy) 


Fig. 17. — Scatter plot of the radio flux densities vs. the y-ray flux for blazar- 
like sources. Spectroscopically confirmed BZBs and BZQs are shown in blue 
circles and red squares, respectively, while ben sources are black diamonds. 
Sources with a radio counterpart in the NVSS are in the upper panel while 
the sources in SUMSS are in the lower panel. 

that these UGSs tend to be brighter than the BZB population 
in the gamma-rays. It is unlikely that this small UGS sub¬ 
sample is constituted by blazars since given their y-ray fiuxes 
Sy and assuming the typical radio-y-ray fiuxes of the blazar 
population, the expected radio flux densities should be greater 
than - 50 mJy for a large fraction, well above the NVSS and 
the SUMSS flux thresholds. On the other hand we also con¬ 
clude that their steep values of y-ray spectral indices are not 
compatible with sources emitting via dark matter annihilation 
(e.g., Belikov et al. 2012; Drlica-Wagner et al. 2014). 
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Radio - y-ray luminosities: Northern sample 
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Radio - Y-ray luminosities: Southern sample 



Fig. 18. — Scatter plot of the radio - y-ray luminosities for blazar-like 
sources. Spectroscopically confirmed BZBs and BZQs are shown in blue 
circles and red squares, respectively, while bens sources are black diamonds. 
Sources with a radio counterpart in the NVSS are in the upper panel while 
the sources in SUMSS are in the lower panel. 


Fig. 19. — Scatter plot of the infrared spectral index evaluated using the 
first two WISE bands vs. the y-ray index for the blazar-like sources in our 
merged list of refined Fermi associations. Spectroscopically confirmed BZBs 
and BZQs are shown in blue circles and red squares, respectively, while ben 
sources are black diamonds. 


Infrared - y-ray fluxes 
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Space Flight Center; the SIMBAD database operated at CDS, 
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Fig. 20. — The scatterplot of the infrared - y-ray fluxes (upper panel) and 
luminosities (lower panel) for the blazar-like sources. Spectroscopically con¬ 
firmed BZBs and BZQs are shown in blue circles and red squares, respec¬ 
tively, while ben sources are black diamonds. 
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All-sky unidentified gamma-ray sources 

(IFGLE and 2FGLE with no y-ray analysis flag) 




Fig. 21.— All-sky distribution of the remaining UGSs with no y-ray analysis 
flags in the reflned association list of the Fermi catalogs. 


Unidentified Gamma-ray Sources 

(no radio matches and no gamma-ray analysis flags) 



Fig. 22.— All-sky distribution of the 40 UGSs that do not have any NVSS 
and/or SUMSS radio source within their positional uncertainty regions at 
95% level of confldence and do not present any y-ray analysis flag. We show 
those lying in the NVSS footprint (i.e., declination greater than -40 °) as 
black circles while the others in the Southern hemisphere are indicated as 
orange circles. 
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Fig. 23.— The distribution of the y-ray spectral index (top) and that of the 
energy flux (bottom) for the 40 UGSs (straight black line) that do not have 
any NVSS and/or SUMSS radio source within their positional uncertainty re¬ 
gions and without y-ray analysis flags in comparison with those of the BZBs 
(dashed blue line) and the BZQs (straight red line). Cumulative distributions 
are also reported on the right panels. The energy flux is the one reported in 
both the IFGL and the 2FGL catalogs. 
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TABLE 7 

Unidentified Gamma-ray Sources with no radio matches and no y-RAY 
ANALYSIS FLAGS. 


IFGL 2FGL 

name name 


1FGLJ0032.7-5519 
1FGLJ0212.3+5319 
1FGLJ0545.6+6022 
IFGL J0709.0-1116 
IFGL J0802.4-5622 
1FGLJ0854.6-4504 
1FGLJ0933.9-6228 
1FGLJ1036.2-6719 
IFGL Ill 17.0-5339 
1FGLJ1306.4-6038 
1FGLJ1405.5-5846 
1FGLJ1441.8-6404 
1FGLJ1518.0-5233 
IFGL J1639.5-5152 
1FGLJ1650.3-5410 
1FGLJ1702.4-5653 
1FGLJ1743.5-3314 
1FGLJ1806.2-H0609 
1FGLJ1825.7-1410C 
1FGLJ1948.6-H2437 
IFGL J2112.5-3044 
1FGLJ2227.4-7804 
1FGLJ2325.8-4043 
1FGLJ2333.0-5535 


2FGLJ0032.7-5521 

2FGLJ0212.1H-5318 

2FGLJ0545.6-H6018 

2FGLJ0802.7-5615 
2FGLJ0854.7-4501 
2FGLJ0934.0-6231 
2FGLJ1036.1-6722 
2FGL Ill 17.2-5341 
2FGLJ 1306.2-6044 


2FGLJ1518.4-5233 


2FGLJ 1702.5-5654 

2FGLJ1805.8-H0612 


2FGLJ2112.5-3042 


2FGLJ2332.5-5535 
2FGLJ0237.9-H5238 
2FGLJ1016.4-4244 
2FGLJ 1058.7-6621 
2FGLJ 1208.5-6240 
2FGLJ 1404.0-5244 
2FGLJ1410.4-H7411 
2FGLJ 1422.3-6841 
2FGLJ 1507.0-6223 
2FGLJ1617.3-5336 
2FGLJ 1626.4-4408 
2FGLJ 1639.8-5145 
2FGLJ1641.8-5319 
2FGLJ1643.3-4928 
2FGLJ1741.1-6750 
2FGLJ 1742.5-3323 
2FGLJ1753.8-4446 


Research in Astronomy and the Leiden Observatory. We ac¬ 
knowledge the WENSS team consisted of Ger de Bruyn, Yuan 
Tang, Roeland Rengelink, George Miley, Huub Rottgering, 
Malcolm Bremer, Martin Bremer, Wim Brouw, Ernst Rai- 
mond and David Fullagar for the extensive work aimed at pro- 
ducing the WENSS catalog. TOPCAT^o (Taylor 2005) for the 
preparation and manipulation of the tabular data and the im¬ 
ages. The Aladin Java applet^ ^ was used to create the finding 
charts reported in this paper (Bonnarell et al. 2000). It can be 
started from the CDS (Strasbourg - France), from the CFA 
(Harvard - USA), from the ADAC (Tokyo - Japan), from the 
lUCAA (Pune - India), from the UKADC (Cambridge - UK), 
or from the CADC (Victoria - Canada). 
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Fig. 24.— Upper panel) The optical spectra of the counterparts associated with 2FGL J1719.3+1744, 2FGL J1801.7+4405 and IFGL J1942.7+1033 already 
available in literature and re-observed at OAN in San Pedro Martir (Mexico) on 29 June and 30 June 2014. The Mgii and the [Oil] 43727 emission lines 
superimposed to the optical continuum are visible in 2FGL J1801.7+4405 and allowed us to classify the source as a quasar, while the lack of features in the other 
two spectra led toward a BL Lac classification. (Lower panel) The 5' x 5' finding chart from the Digitized Sky Survey (red filter). The potential counterparts 
pointed during our observations is indicated by the red circle. 


APPENDIX 

Here we present the spectra of 2FGL J1719.3+1744 and 2FGL J1801.7+4405, both with uncertain redshift, as reported in 
the Roma-BZCAT together with that of IFGL J1942.7+1033. We confirmed the redshift for 2FGL J1801.7+4405 and the 
BL Lac classification of IFGL J1942.7+1033 (Tsarevsky et al. 2005; Masetti et al. 2013) while 2FGL J1719.3+1744 (alias 
BZBJ1719+1745) has a completely featureless spectrum. 



























